Fluctuations of radiation from a chaotic laser below threshold 
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Radiation from a chaotic cavity filled with gain medium is considered. A set of coupled equations 
describing the photon density and the population of gain medium is proposed and solved. The 
spectral distribution and fluctuations of the radiation are found. The full noise is a result of a 
competition between positive correlations of photons with equal frequencies (due to stimulated 
emission and chaotic scattering) which increase fluctuations, and a suppression due to interaction 
with a gain medium which leads to negative correlations between photons. The latter effect is 
responsible for a pronounced suppression of the photonic noise as compared to the linear theory 
predictions. 
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Introduction. Amplifying random media have now 
been a focus of intense research for almost a decade 
[1 H H II 13 The ultimate goal of these efforts is a 
construction of a random laser in which laser action uti- 
lizes the feedback provided by multiple scattering from 
disorder. Lasing action was reported in a recent work by 
Cao et al. 0, Q and Frolov et al. The difference 

between the true coherent random laser (RL) and the 
amplified spontaneous emission (ASE) cannot be identi- 
fied by the emitted light intensity or its frequency profile 
as both ASE and RL are characterized by a pronounced 
narrowing of the spectral linewidth 7 compared to the 
atomic transition width T. The knowledge of the statis- 
tics of light fluctuations is therefore necessary to make 
a distinction. While RL action should be coherent and 
therefore obey the Poissonian distribution, the ASE fluc- 
tuations are above the Poissonian level (super-Poissonian 
noise). Both cases seem to be realized in experiment, that 
is of RL and ASE [13. There is still no complete 

theory able to predict which case ought to be expected 
in a particular system. From the outset, one expects RL 
regime to occur whenever there is a single mode above 
the laser threshold or a few non-overlapping modes, and 
ASE case to be realized as long as there is a number of 
overlapping modes above the threshold. In the latter case 
the scattering between the lasing modes should lead to 
the enhancement of fluctuations. Since the broadening 
of the laser modes is given roughly by the inverse mean 



dwell time 



dw ' 



cavity modes overlap strongly provided 



that the mean level spacing A is small, rlj 3> A. (Here- 
inafter, we use the units with h = 1.) 

Although there is by now a substantial literature on 
various models for random lasers 0, I12I IT^L IlI ITfij , we 
address only those papers that consider fluctuations of 
laser radiation. Beenakker put forward a quantum ap- 
proach for a photon statistics based on the random ma- 
trix theory [16j for a linear amplifier below laser thresh- 
old. The threshold is achieved when the dwell time inside 
the system t^w exceeds the characteristic amplification 
time determining the rate for photon emission r a . The 
threshold features a singularity in the photon density. 
The fluctuations of radiation demonstrate even stronger 



singularity as the ratio of the noise power to the mean 
photon flux S/J diverge at the threshold, as a result of 
scattering between many lasing modes. Linear theory 
fails upon approaching the threshold when the assump- 
tion of equilibrium gain breaks down. Hackenbroich et 
al. considered statistics of a single lasing mode in a 
chaotic cavity above the threshold. Cavity opening was 
assumed to be covered with a barely transparent mir- 
ror making the inverse dwell time small compared to the 
spacing between the modes, <C A, thus preventing 
them from overlapping. In Ref. [Tsj the nonlinear ef- 
fects have been numerically studied for a system with an 
opening much smaller than the wavelength. While such 
a model might be useful for a microwave emission, its 
relevance to random lasers with typical openings much 
wider than the wavelength of light remains unclear. 




FIG. 1: Chaotic laser: a cavity characterized by the mean 
level spacing A contains A^o particles of a gain medium and 
is connected to the outside via a waveguide (lead) supporting 
M transverse channels. 

Chaotic Laser. In the present paper we develop a semi- 
classical approach based on the radiative transfer theory 
|l9j | for a chaotic cavity filled with a gain medium, see 
Fig. 1 (let us, following Refs. 01 an d [13 ■ call such a 
system 'chaotic laser'). The cavity is coupled to the out- 
side through a waveguide supporting M 3> 1 transverse 
modes, so the inverse dwell time tT 1 — -^A is large com- 

7 aw ztt 

which is in- 



pared to the mean level spacing A 

versely proportional to the photon density of states and 
a system volume L 3 . The ergodic time T er = L/c is the 
shortest time scale of the system, r er <C Tdw , T a . This 
condition assures that photons scatter chaotically from 
the boundaries many times before leaving the cavity or 
stimulating an emission of another photon, therefore re- 
sulting in the homogeneous and isotropic distribution f u 



2 



everywhere inside the cavity. The number of photons in- 
side the cavity with the frequency in the interval dtu is 
given by f^doj/A. The rate at which this distribution 
changes in time, 



at 



(1) 



is determined by the photon flux via the waveguide, J u 
and by the photon creation and absorption inside the 
cavity I u . The flow of photons through the waveguide 
(lead) is proportional to the size of the opening (charac- 
terized by the number of channels M) and the difference 
between the concentration of photons inside and outside 
the cavity, the latter denoted by n u (in general case of 
external radiation incident on the system), 



Z7T 



(2) 



with the last term standing for the Langevin source de- 
scribing fluctuations in the system. The left-hand side of 
the rate equation Q has a familiar form of the particle 
conservation condition while the right-hand side accounts 
for the photon creation and absorption inside the cavity, 



I u = wZ(fuj + 1) - w u f u + C W A. 



(3) 



The quantities u>+ and w~ stand for the rates of photon 
emission and absorption respectively. The probability 
of photon creation w^ J (f LJ + 1) takes into account both 
stimulated and spontaneous emission as required by the 
quantum statistics of photons. 

In addition, Eqs. lll.'il) contain stochastic Langevin 
sources which are necessary to take into account the ran- 
domness of photon transmission through the lead J u , and 
the randomness of the emission and absorption events, 
£ u . These terms have zero average, J u = £ u = 0, 
and correlators which follow from the assumption that 
elementary stochastic events of scattering, emission and 
absorption have independent (and therefore Poissonian) 
distributions. In particular, the Langevin term associ- 
ated with non-conservation of particles, C u = + £j 
consists of separate contributions from emission and ab- 
sorption. Different processes have independent distribu- 
tions (in particular, they are uncorrelated when occur at 
different times) with the second moments given simply 
by their average rates (property characteristic for Pois- 
sonian processes), 



£S(*)C(*0 = tUA'^+^ + i), 

£J (*)£", (tf) = S^tA-'wZJu, £+(*)£-,(*') =0, 

where <5 wt is the shorthand for S(ui — u')8(t — t'). Using 
these formulas we arrive to the following expression, 



CuW u .(t') = <UA->+(/ w + i) + w -f u 



(4) 



Similar reasoning allows us to write the correlation func- 
tion for the Langevin sources describing randomness in 



the photon flux through the waveguide, 



A M \~ (- 
Ztt 



+ + (5) 



Since we assume no reflection in the lead, the only re- 
maining source of fluctuations is the randomness of the 
photon densities at both ends of the waveguide [2(j . 

In equilibrium at temperature T, the function f u must 
be equal to the Bose-Einstein distribution. This fixes 
the ratio, w~/w+ = e"/ T . So far our discussion has 
been completely general, valid for both the absorbing, 
W Z > w t i T > °)> and amplifying < w+ (T < 0), 
media. Hereinafter we concentrate for simplicity on the 
case of complete population inversion, w~ = (T = 
0~). The complete population inversion is realized, e.g. 
when the radiative transition occurs not to a ground state 
but rather to some excited state with very short (non- 
radiative) lifetime, see Fig. 2. 



pumping 



radiative transition 



fast non-radiative 
transition 







FIG. 2: Complete population inversion realized in a three- 
level laser. The atom of the gain medium is pumped to an 
excited state 2, the photon is emitted when the atom under- 
goes a transition to a lower excited state 1. If the relaxation 
from the excited state 1 to the ground state is fast, the ex- 
cited state 1 is almost always empty. Therefore the radiative 
transition 2— >1 occurs with the effective temperature T = CP. 

We assume the spectral line for atomic transitions to 
be a Lorentzian centered around the frequency ujq with 
the linewidth T, 



;+ — 



jNA, 



WqT 



n[{u - cj ) 2 + T 2 ] 



(6) 



where N is the number of excited atoms of the gain 
medium in the cavity, and Wq is the total probability (per 
unit time) for an excited atom to spontaneously emit a 
photon into vacuum. 

The dynamics of the gain medium follows the equation, 



Ot T 



dujI^+Af, (7) 



where P is the pumping rate, Nq is the total number of 
gain medium particles, and r _1 is the probability for a 
non-radiative relaxation. The third term is the total loss 
from radiative transitions into all possible cavity modes. 
Note that it contains stochastic emission source, as spec- 
ified by Eq. ©. The Langevin term N describes the 
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spontaneous fluctuations in the number of excited atoms, 
which are also assumed Poissonian, 



N{t)N{t') = S(t - t')[P(N -N) + N/t] 



(8) 



Note that different Langevin terms C u , Af are in- 
dependent and have zero cross-correlators. The coupled 
equations , together with the correlators 1 1131) , © 
determine the radiation from a chaotic cavity. 

Mean photon distribution. We illustrate the above for- 
malism by calculating the average photon flux from the 
cavity and its fluctuations assuming no radiation inci- 
dent on the system, n u — 0. In the stationary regime the 
time-averaged equation for the gain medium Q is easily 
solved and the number of excited atoms is expressed via 
the photon distribution, 



N 



PN Q 



Pi + w + J dw a u f u 



(9) 



where P\ = P + r — 1 . Substituting it into Eqs. i|H3|) one 
obtains an integral equation for the photon distribution 
function, which could be solved analitically. The time- 
averaged distribution has the form of a Lorentzian, 



f LJ 



r 2 - 7 2 



(cj - loq) 2 + i 2 



(10) 



with the lincwidth 7 determined from the algebraic equa- 
tion, 



r 2 - 7 2 



2N w TjP 
M(w T + P17) ' 



(11) 



The total photon flux leaving the cavity J = || J duof u 
is found from Eqs. I|10ll lfl . 



J = 



NqWqTP 
wqF + Pi 7 



(12) 



The cubic equation Ijlll) allows easy numerical solution 
for any values of the pumping strength P and the num- 
ber of channels M. Here we analyze in detail the more 
illuminating case of absent non-radiative losses, r _1 = 0. 

It is convenient to introduce the dimensionless pump- 
ing strength and the effective number of radiative tran- 
sitions, 



P 

P = — , 
wo 



v 



2N w 

my 



The parameter rj is the ratio of the photon dwell time in 
the cavity, Td w = 2ir/MA, and the characteristic time of 
spontaneous photon emission, which at the center of the 
line is given, according to Eq. ©, by r a = nT/N woA 
(if all excited states are filled). The dimensionless width 
g = 7/r satisfies a simple equation, 



1 



.9 



tjpg 



P9 



(13) 



This equation determines the radiation linewidth g v (p) 
as a function of the pumping strength for different sys- 
tem geometries: with the lower values of rj characterizing 
open cavities (large number of outgoing channels M) and 
the higher values of r\ corresponding to pinching off the 
waveguide opening. The other two relevant quantities, 
the population of the gain medium and the outgoing flux 
can also be expressed in terms of the new dimensionless 
variables, 



N 



PfJ 



2J_ 

MT 



VP 



(14) 



1+pg MT 1 + pg 

For weak pumping, pC 1, the width is narrowed linearly 
with the pumping strength, 

VP 
2 ' 



1 



(15) 



with the narrowing pronounced stronger for less open 
cavities (larger values of 77). The strong pumping, p> 1, 
behavior is qualitatively different for different values of 
the effective number of transitions rj. 

(i) Open system, ?y < 1 (large number of channels sup- 
ported by the waveguide). The linewidth remains finite 
for any value of the pumping strength, with, 



lim g v (p) = y/T~^ 



V- 



The population of excited atoms becomes saturated, 
N — > No, while the outgoing flux reaches a constant 
value J oc (1 — 77) _1 ^ 2 . The flux diverges upon pinch- 
ing off the waveguide opening and approaching the point 
rj = 1. This point is referred to as the 'laser threshold' 
in the linear amplifier theory [l6l |. However, it is not a 
real laser threshold (in a sense discussed towards the end 
of the paper) but rather a point where nonlinear effects 
can no longer be disregarded. We call it a 'critical point' 
here. 



(ii) In the critical point, rj = 1, the gain population 
remains full at large pumping rate and the outgoing flux 
diverges: J oc p 1 / 3 . The linewidth is a slowly decreasing 
function of pumping g oc p" 1 / 3 . 

(iii) Closed system, 77 > 1 . In this case the linewidth 
decreases inversely proportional to the pumping inten- 
sity. The number of photons leaving the cavity grows 
linearly with the pumping. The gain population remains 
partially empty: N/N — > 1/rj < 1. The results for the 
strong pumping regime are summarized in Table 1. The 
quantitative behavior of the linewidth is illustrated in 
Fig. 3. 

Fluctuations. We now consider the fluctuations of ra- 
diation around the mean solution obtained above. Both 
the photon distribution and the number of excited atoms 
deviate from their mean values, f u = f u + 6f u , N = 
N + 5N. The fluctuation of the number of excited atoms 
is found from the linearized equation Q, 



SN = N 



M — J duj£ u — N j duo a^Sfu, 
PiVo 



(16) 



TABLE I: The dependence of the dimensionless linewidth g = 
7/T, population of the gain medium, and the outgoing photon 
flux for large values of the pumping intensity p = P/wo ^ 1, 
for different values of the effective number of optic transitions 
77 (in the absense of non-radiative losses). 
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77 = 1 


1 

pT/3 
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1/3 


77 > 1 


1 

(7? - l)p 


1 
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Substituting it into the equation we obtain the linear 
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FIG. 3: The linewidth of radiation from a chaotic laser, 
g — 7/r, versus the dimensionless pumping intensity p = 
P/wq (neglecting non-radiative losses, t" 1 = 0) for differ- 
ent values of the effective number of radiative transitions, 
77 = 2Now /Mr. Solid line: 77 = 1/2, the linewidth g tends 
to the limit l/y/2 w 0.71 for p — > 00. Dashed line: 77 = 1, 
the linewidth decreases as g w p~ 1,/3 . Dotted line: 77 = 2, the 
linewidth behaves as g w p _1 . 



integral equation for the fluctuation of the distribution 
function of photons, 



Sf u 



Mil 



A/" - / £ a 



2ira u PN J ^ JW ttuJ N PN a 

Solving this integral equation is rather straightforward and yields the fluctuation of the outgoing flux, 



6J u = ^5f u + J u = (f u + - f u J u - A, I ,U (f _ y -- I /.,X 



(17) 



(18) 



In the last expression we utilized the following notations, 



2nPN D ' 271-PiVo 



By using the correlation functions of the Langevin 
sources Q-© an d iJHJ the frequency-resolved noise 
power, S UU ' — J dt J u (t)Ju'(0), is found, 

Su)uj' = S(u - ui')S u — {A u Su' + A'Suj) 

+A^A lj ,[S + P{N -N)+t- 1 N}, (20) 

where we introduced the notations, 



S u = JM„ + !)(/« + fu, + 1), S = J duS u . (21) 

The first term in Eq. H20|) represents the correlations be- 
tween the photons of the same frequency. Those cor- 
relations are the result of spontaneous and stimulated 



emission and also due to chaotic photon scattering inside 
the cavity. Assuming again the absence of non-radiative 
losses we can express the Fano factor (the ratio of the 
noise power to the mean flux) at the center of the line, 
u> = u)q, via a single parameter, 



TP _ j^O _ 

*o-=-- 

• J CJn 



1-ff 2 



9 



a" 



(22) 



namely the dimensionless linewidth g — 7/r. This 
is a monotonously increasing function of the pumping 
strength p. As before, we analyze three possible regimes, 
making use of the above results summarized in Table 1. 

(i) Open system, 77 < 1. When the pumping intensity 
increases, p — > 00, the Fano factor approaches the finite 
value, 



F n = 



1 — 77 + 77 

(1 — v) 3 



(23) 



This expression reproduces the result of the random ma- 
trix theory |16( for a linear amplifier with a complete 
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population inversion. Indeed, in an open system, r\ < 1, 
the gain medium dynamics becomes suppressed for large 
values of the pumping strength when all excited states 
become occupied. 

(ii) Critical point, rj = 1. For large pumping intensity, 
p ^> 1, the Fano factor grows quadratically, 

F =P 2 . 

(i) Finally, in a closed system, rj > 1, the Fano factor 
behaves as the sixth power of the pumping strength p for 
P> 1, 

F = (r, - 1) P 6 . 

The quantitative behavior of the Fano factor (12211 is il- 
lustrated in Fig. 4 where numerical results are presented. 

Cross-correlations. Contrary to the linear theory pre- 
dicting the photons of different frequencies to be uncor- 
related, we observe significant off-diagonal correlations, 
as attributed by the expression l|20|l . This is due to the 
strong coupling of photons to the fluctuations of the gain 
medium. The important fact is that these correlations 
are negative. It could be understood from the following 
reasoning: the emission of a photon at some frequency 
reduces the population of excited atoms making therefore 
the emission of another photon less likely. 

To illustrate better the role of cross-correlations it is 
useful to calculate the noise power for the total radiation 
from the cavity. Integrating the expression l|2U[l over both 
frequencies yields the noise power S — J duduj' S^^i for 
the total photon flux leaving the system, 

S = JA 2 + +S{l-Af, (24) 

\ w Tt J 

here A — J dioA^. By using the expressions (|§|)- (|10fl for 

the mean distributions, the necessary integrals are easily 
evaluated, 

~ = M(T* - 7 2 ) (5r6 _ 3F 4 72 + + ? 6) 

327' 

_ M(r 4 - 7 4 ) M(r 2 - 7 2 ) 2 

APNq^D ' 4PiV 7 3 ' 

Again, we analyze the obtained general expression ne- 
glecting non-radiative losses, r _1 = 0. By using the di- 
mensionless variables p and 77 defined in the preceding 
sections, the Fano factor for the total photon flux leav- 
ing the cavity, F — S/ J, can be written as, 

4(l + g 2 ) 2 + P 2 (5-3 g 2 + 7, 9 4 + 7 g 6 ) 

4(1 + <? 2 + 2p 5 3 ) 2 ' 1 J 

The dimensionless width g is a function of both the 
pumping strength p and the number of optic transitions rj 
and is given by the solution of the equation l|13l) . The de- 
pendence of the Fano factor (|25J) on the pumping strength 
is plotted in Fig. 5. A surprising feature of this depen- 
dence is the suppression of the noise below Poissonian 



level for intermediate values of the pumping strength. 
This is the demonstration of negative correlations be- 
tween the photons due to their interaction with the gain 
media. Indeed, despite the fact that the diagonal corre- 
lations (|21I22I) increase monotonously with the pumping 
p as a result of photon bunching, the total flow of out- 
going photons becomes less fluctuating. For large val- 
ues of p the enhancement of the photon bunching over- 
whelms the negative correlations and the Fano factor F 
increases. Still, it remains strongly reduced compared to 
the frequency-resolved Fano factor (|22|l . see Fig. 5. 

(i) For an open system, rj < 1, and large values of the 
pumping strength the Fano factor approaches the limit, 

1-2^+^-^3 
(I-,?) 3 

and diverges at the critical point rj — > 1. Its ratio, F/Fo, 
to the frequency-resolved Fano factor (|23[) remains con- 
stant and decays monotonously from 1 to 5/16 while rj 
changes from to 1. 

(ii) Critical point, r\ = 1. The Fano factor increases 
quadratically for large pumping, p> 1, 

36 ' 

with the ratio of the two factors remaining finite, F/Fq — » 
5/36. 

(iii) Finally, in a closed system, 77 > 1, the fluctuations 
of the total flux are strongly suppressed, 

4 ' 

compared to the diagonal fluctuations, with the ratio 
F/F a -» 0. 

Laser threshold. So far the mean level spacing A has 
not entered the expressions for the mean flux nor the 
noise power. However, it sets the condition for applica- 
bility of the entire approach. The above equations are 
semiclassical and imply incoherent radiation from the 
system. This is so provided that there are many over- 
lapping modes inside the cavity within the linewidth 7. 
A characteristic mode width is given by the inverse dwell 
time r^J = 4^-. As long as 7, ttJ 3> A, the mode over- 
lapping is strong. Estimating the linewidth as 7 ~ T/r/p 
(for rj > 1), we can establish the necessary conditions: 
» A, and M > 1. 

Another condition for the semiclassical theory to be 
valid requires the system to remain below the laser 
threshold. The latter condition, for the i-th. cavity mode, 
is determined by the relation, Tj^/to = = 1, with 

(i) 

the amplification time Ta being in general a function of 
the photon population of all cavity modes and, hence, 
some complicated function of the pumping strength. 
From our equations © and J§J) the semiclassical value 
of s at the center of the line uj = luq is found, 

s = -^ = l-g 2 <l, 
a^NA 
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and the system stays below the threshold for any pump- 
ing intensity. However, as the pumping increases, g — > 0, 
the quantum distribution of around the mean value 
s becomes important |2l| . At some critical pumping, the 
mode with the largest s^' reaches its threshold value and 
starts to lase in a coherent fashion. This regime is be- 
yond the reach of the semiclassical theory. The problem 
of fluctuations of radiation from a chaotic laser above the 
threshold in the non-linear quantum regime remains an 
open one. 

Conclusions. We presented a semiclassical theory of 
the incoherent radiation from a chaotic cavity containing 
optically active medium based on the radiative transfer 
theory with Langevin sources. The theory reproduces the 
results of the random matrix approach for fluctuations of 



radiation from a linear amplifier. By accounting for the 
dynamics of the gain media the theory is extended to in- 
corporate the nonlinear effects: in particular, the spectral 
line narrowing of the emerging light is calculated. The 
fluctuations of radiation from the system are the result of 
two competing effects: the enhancement of fluctuations 
due to photon bunching and the suppression of fluctua- 
tions due to interaction with the gain medium. 
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FIG. 4: The Fano factor for the fluctuations of radiation at the 
center of the line, u) — ljq, versus the dimensionless pumping 
intensity (neglecting non-radiative losses) for different values 
of the effective number of radiative transitions: r\ — 1/2 (solid 
line), 1 (dashed line), 2 (dotted line). 



FIG. 5: The Fano facor for the total radiation from a chaotic 
laser for different values of the effective number of radiative 
transitions: ij = 1/2 (solid line), 1 (dashed line), 2 (dotted 
line) . 



